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1 Introduction 

Since Kowalevskaya's monumental work [15], the Painleve test has been the 
most widely used and the most successful technique for detecting integrable 



differential equations [H IH E]- The test has been apphed to many differ- 
ential equations and, for those passing the Painleve test, the indicators of 
integrability such as the existence of sufficiently many conservation laws, the 
Lax pair, the Darboux transform, the Backlund transform have always been 
found. 

The Painleve test itself is a recipe for finding all formal Laurent series so- 
lutions with movable singularities. An n-th order system passes the Painleve 
test if such formal Laurent series solutions admit n free parameters (including 
the location to of the the movable singularity as one free parameter). Since 
the solutions of an n-th order system should have n degrees of freedom, the 
Laurent series solutions, with maximal number of free parameters, should 
include all the solutions. Then one hopes that all movable singular solu- 
tions are poles and concludes that all solutions are single valued. This is the 
heuristic reason why passing the Painleve test means integrability. Indeed 
such reason underlies the classification of integrable equations by Painleve 
and the others [T6] . 

In [21 E], Adler and van Moerbeke put the heuristic reason on solid foun- 
dation for the very nice case of algebraically completely integrable systems. 
They used the toric geometry to construct a complete phase space, and 
thereby gave a satisfactory explanation for the relation between the Painleve 
test and the integrability. In [5], Ercolani and Siggia pointed out that al- 
gebraic geometry is not needed for integrable Hamiltonian systems. They 
suggested using the expansion of the Hamilton- Jacobi equation to construct 
the change of variable used for completing the phase space. They showed 
that the method often works through many examples, but did not prove that 
the method always works. 

In this paper, we show that, not only is algebraic geometry not needed, 
the Hamiltonian set up by Ercolani and Siggia is also not needed. In fact, 
with virtually no condition, the existence of the change of variable needed 
for completing the phase space is equivalent to passing the Painleve test. 

Theorem A. A regular system of ordinary differential equations passes the 
Painleve test if and only if there is a triangular change of variable, such that 
the system is converted to another regular system, and the Laurent series so- 
lutions produced by the Painleve test are converted to power series solutions. 

For the exact meaning of "passing the Painleve test", see the rather 
straightforward definition in the beginning of Section |2l Also see jTT] for 



some concrete examples. The condition is as straightforward and elemen- 
tary as can be. Although the proof is also rather elementary, we believe 
there is some interesting underlying algebraic structure that is worth further 
exploration. 

The theorem gives a key connection between the two major works by 
Kowalevskaya [13]. As in the theorem of Cauchy and Kowalevskaya, we 
consider a regular system of differential equations with complex analytic 
functions on the right side. Given a movable singular solution, represented 
by Laurent series in (t — to) for variables Ui, changing Ui to u~^ certainly 
regularizes the variables. However, such simple regularization will create 
singular differential equations for the new variables. By the resolution of 
singularity, we mean a change of variable that regularizes both the solution 
and the equation. This is not always possible, and the condition for the 
resolvability of the singularity is exactly passing the Painleve test. 

The regularization is achieved by a triangular change of variable 

ui = T-\ 

Ui = ai{t, T,p2,..., Pi-i) + bi{t, r, p2, • • • , Pi-i)pi, I <i <n, 

where —k is the leading order of the Laurent series solution for ui obtained 
in the Painleve test, ai,bi are meromorphic in r and analytic in the other 
variables, and bi does not take zero value. Such change of variable can be 
easily inverted, at the cost of introducing one k-th root. Then it is easy to 
see that the Laurent series for ui, . . . ,Un correspond to a power series for 
r satisfying r(to) = 0, r'ito) ^ 0, and the Laurent series for p2, • • • ,Pn- In 
case the system for ^i, . . . ,«„ passes the Painleve test, we can find suitable 
ai,6j so that the Laurent series for p2^...^Pn are always power series. In 
fact, in our change of variable formula f l2.3p and fl2.18p (or f l2.3p and fl2.12p ). 
we get very specific forms for Oj and b^. Then applying Cauchy theorem to 
the regular system for r, /92, • • • ,Pn and the convert to the original system, 
we conclude the following. See [TJ Section 3.8.6] for more discussion on the 
convergence of formal Laurent series solutions. 

Corollary. // a regular system of ordinary differential equations passes the 
Painleve test, then the Laurent series solutions obtained from the test always 
converge. 

Section |2]is devoted to the rigorous proof of Theorem A. We will be rather 
conservative in assuming that the right side of the differential equations are 



polynomials. The assumption is used to make sure that, after substituting in 
the Laurent series, the right side again becomes a Laurent series. Certainly 
we can make sense of this with functions other than polynomials, but the 
spirit of Theorem A is already fully reflected in the polynomial case. 

Ercolani and Siggia's work suggests that there should be a version of 
Theorem A that is compatible with Hamiltonian structure. 

Theorem B. // a regular Hamiltonian system of ordinary differential equa- 
tions passes the Painleve test in the Hamiltonian way, then there is a canoni- 
cal triangular change of variable, such that the system is converted to another 
regular Hamiltonian system, and the Laurent series solutions are converted 
to power series solutions. Moreover, if the system is autonomous, then the 
new Hamiltonian function is obtained by substituting the new variables. If 
the system is not autonomous, then the new Hamiltonian function is obtained 
by substituting the new variables and then dropping the singular terms. 

The exact meaning of the "Hamiltonian way" is given in Section [31 See 
fr2\ for some concrete examples. 

In view of the rather ad hoc nature of the Painleve test, we further provide 
a rigorous foundation of the Painleve test in Section |H Our setup is by no 
means the broadest possible (and is therefore only "a" foundation), but is 
satisfied by most examples and is the most common way the conditions of 
our theorems are met. 

We emphasize that the results of this paper are local, in fact as local as the 
theorem of Cauchy and Kowalevskaya. Although the Painleve test is local, it 
is supposed to test the global properties of solutions. The regularization in 
this paper may be used to further prove the global property for some specific 
integrable equations, but here we are not making any general statement about 
the global property. See [HI El [101 [IH HZ] for specific examples of global 
Painleve analysis. 

Finally, we remark that our theorems do not deal with the lower balances 
between the principal balances, as discussed in [21 [3l [5] as part of the process 
of completing the phase space. If there are no lower balances, then our the- 
orems effectively show that passing the Painleve test implies the completion 
of the phase space and thereby explains the integrability. Further extensions 
of our theorems to lower balances are needed for the more general case. 



2 Resolution of Movable Singularity 

Consider a system of ordinary differential equations 

fi{t,Ui,...,Un), 



U, 



1 < i < n. 



We restrict the discussion to fi(t,ui, . . . ,Un) being analytic in t and poly- 
nomial in ui, . . . ,Un- At the end of this section, we will point out what is 
exactly needed for more general /j. 

A balance for the system is a formal Laurent series solution, with movable 
singularity to- The balance often admits several free parameters, similar to 
the initial values in formal power series solutions. Therefore a balance is 
typically given by (ri is reserved for the first free parameter to) 



Ui = a,,o(to)(t - to) '^^ + ■ ■ ■ + ai,Ai^i(to)(t - to)^^ 



Ai-l-fci 



+ ai,Ai(to,r2,...,r„J(t-to)^^ ''^ 
H \-ai^x2-i{tQ,r2,...,rn^){t 

+ ■■■ 



\A2-1-fei 



(2.1) 



+ ai,A.(to, ^2, . . . , r„J(t - to)-^^' 
H haij(to,r2,...,r„J(t 



toY-"^ 



+ 



where aij is analytic in its variables. The free parameters are called reso- 
nance parameters, and the indices Ai, . . . , A^ where they first appear are the 
resonances. In the balance (12. ip . the resonance parameters r„j_j+i, . . . ,rni 
have the corresponding resonance A;. The number nii = ni—ni_i of resonance 
parameters of resonance A; is the multiplicity of A;. We have 



ni 



mo + mi + \-mi, 



where mg = 1 counts the resonance parameter to of resonance Aq = —1. 

We certainly want the free parameters to really represent the variety 
of Laurent series solutions. This means that we should require that the 
resonance vectors of resonance A/ to form an n x tti; matrix of full rank nii 



/ da 



Ri 



1,A, 



dr 



dttn^Xi 



dai,xi \ 
dr. 



ni 



da. 



n,X, 



\d 



ni-1+1 



drn, J 



For the resonance Aq = —1, this means that the basic resonance vector 
Rq = {—kiOifi, . . . , —knanfl)^ is non-zero. In fact, when we consider the 
independence among all the resonance parameters, we really should require 
that the resonance matrix 



R = [Rq Ri 



-Rs 



to have the full rank Ug- 

The exact condition for Theorem A is the following. 

Definition. A balance (12. ip is principal, ii rig = n and the resonance matrix 
is invertible. 



2.1 Indicial Normalization 

By substituting to = ^ — (^ — ^o), the coefficients Uij become power series 
of {t — to) with analytic functions of t,r2,...,r„ as coefficients. Further 
substituting these power series into (12. ip . the balance becomes 



.(0)/ 



(0) 



u, = a^>it)it - 1,)-'^ + ...+ a^i^_,mt - tor^ 



\\]_-l-ki 



SO) 



Xl-ki 



+ ■ ■ ■ + «S.-i(^' ^2, . . . , r„J(t - to)"^-^-'=" 

+ ■■■ 

+ aS(t,r2,...,r„J(t-to)^=-'^' + ---. 



We have 



«i,Ai(^''^2, . . . ,r„J = ai,xi{t,r2, . . . , r„ J + Q,A,(t,r2, 
Therefore the new resonance matrix 



i?(°) = (i?r M°^ ■■■^i°^ R. 



(0) 



dr. 



"i-i+i 






\^^m-i- 



'"'"l-iJ- 



dr, 



da. 



ni 



(0) 
n,A, 



d^m j 



(2.2) 



is obtained from the old one simply by replacing tg with t 

The matrix is still invertible. 

Up to rearranging the order of Ui, we may assume kia\ o(t) 7^ 0. Then we 
introduce a new variable r by 

Ui = T-''K (2.3) 

Fixing some {—ki)-th. root /3i(t) of a-j^Q(t) and taking the (— A;i)-th root of 
the equality (12. 2 p for ui, we get the power r-series 

r = /3i(t)(t - to) + ■ • ■ + /3Ai-i(t)(t - toy^-' 

+ PxAt, r2,..., r„J(t - to)^' + ■■■ + I3x,-i{t, rs, . . . , r„J(t - to)^^"^ 
+ ■ ■ ■ (2.4) 

+ /3A.(t,r2,...,r„J(t-to)^= + ---, 
We may reverse the power r-series and get 

t-to = hT + hr^ + hr^ + ■■■ , (2.5) 

where bi is an analytic function of the same variables as /3j. 

Substituting the power r-series (12.51) into the balance (12.21) . we get Lau- 
rent r-series 

u, = af^b^'^^T-^^ + [-hafj,b:,% + af}bi]b-^''^ r^^''^ + ■■■ (2.6) 

+ [-kiaf)^b^%+i + afjh'^ + terms involving aJ^J^-, b<j\b^^'T^-^' + ■■■ . 

For Ml = r^'^i, the equality (12. 6p gives the recursive relation that computes 
the coefficients bj 

bj+i = ^ to) (^i,j + terms involving a[°^j-, b<j. 

"'1^1,0 

Substituting the coefficients bj into the other Ui in (12. 6p . we get updated 
Laurent r-series 

+ a^lit, r2, . . . , r.Jr^-'^^ + ■ ■ ■ + a« _,(t, r^, . . . , r^Jr^-^-^^^ 

+ ■■■ (2.7) 

+ a^Ut, r2, . . . , r^Jr^^-'^^ + ■ • • , z > tiq = 1, 



where 






"S- 



, (0) 



terms involving a\ 



(0) ^(0) 



<ji "'i,<j^ "<j 



(1) ;. 



,&<,-, i > 1. 



For j = Xi, i > 1 and ?7,/_i < s < ni, the formula for a^ implies 



da) 



(1) 



J,A, 



5r^ 






Therefore the updated resonance matrix 



dr 



i?« = (i?l'^ ■■■ i?«), i?! 



(1) 



ni-i+l 






da. 



(1) 

n,A; 



5a! 



(1) 

n,A, 



\dr 



"i-i+i 



^^n. / 



is obtained by deleting the first column (which consists of entirely 0) of the 
following matrix 






b^'-J 



/ ^202,0 

kia^ 



k 0^°^ 

\ M«1,0 



\ 



i?(°) 



"1 J-mi 



V 



^i inis/ 



J 



In other words, up to multiplying the powers of hi to rows and columns, i?'-^-* 
is obtained as part of the row operation that uses the first entry of R^^^ to 
eliminate the other terms in the first column. In particular, the invertibility 
of R^'^'> implies the invertibility of R^^^ . The size is reduced by 1 because the 
resonance parameter tg has been "absorbed" into the new variable r. 



2.2 Resonance Variable 

Next we introduce new variables to "absorb" the resonance parameters r(i „j] = 
(r2, . . . , r„J of resonance Ai. By rearranging U2, ■ ■ ■ ,Un if necessary, we may 



assume that the first mi x mi submatrix A^^'^ of i?'-^^ is invertible. Then we 
have 



where 



^S^^ - (^5(1) j ' (^?^ ■ ■ ■ ^^^^) - (^/}(i; 



We introduce new variables p2, . . . , Pm by truncating the r-series of W2, • • • , li^^ 
in ([22D at r^i"^* 

«. = «S Wr-'^" + ■ ■ ■ + o^Li(t)r''-'-'^ (2.8) 

Then we have the equahties 

«1ai(^'P2, • • ■,pnj = a^ll{t,r2, . . . ,r„J + ■ ■ ■ + ag^„i(t,r2, . . . , r^Jr^^-^-^^ 

+ ■■• 

+ ag^(t, r2, . . . , r„jT^=-^^ + ■ ■ ■ , 1< z < m. 

By the inverse function theorem, the invertibihty of 



A«(t,r2,...,r„J 



5('^2,---,r„J 



means that the left side is a locally invertible map of P2, ■ ■ ■ ,Pnu and for 
small r, the right side is a locally invertible map of r2, . . . , r„^. By taking the 
inverse of the left side map, we get the power r-series for the new variables 

Pi = ri + f3i^i(t, r2, . . . , rnjr H 

+ A,A2-Ai (t, r2, . . . , Tni, r„i+i, . . . , r^Jr^'^^' H 

+ ... 

+ A,A.-Ai(t,r2,...,rn,,rn,+i,...,rnJr^''"^^ H , 1 < i < ni. 



By taking the inverse of the right side map, we get the power r-series for the 
resonance parameters 

Ti = Pi + bi^i{t, p2, . . . , pm)r ^ 



+ 6i,A2-Ai(t,p2,...,Pni,rni+i,...,r„Jr^2 ^1 + .. 

+ ■■■ 



(2.9) 



+ &i,A,-Ai(t,P2,---,Pni,r„i+i,...,r„jT = ^H , l<i<ni. 



Substituting (12. 9p into the balance (12.71) . we get the Laurent r-series 
u^ = aS(t)r-^^ + ■ ■ ■ + a^l_,{ty^-'-'^ + a^l{p,, . . . , p^Jr^-'^' 



+ 

+ ■■ 
+ 



da':'^ 



i,Ai 



(1) 



d{r2,...,rn,) 



{t,p2,..., pm)hi + a^x^+i{t, p2, • • • , Pm) 



„Ai + l-fci 



(2.10) 



"^iM /^ _ _ M , (1) 



^, — At, P2, ■ ■ ■ , Pnr)hj + alx^+j{t, P2, • • • , Pm, r^i+i, • • • ) 

0\T2, • • • ) ^ni j 

+ terms involving a\ <Ai+i' ^<j 



^Ai+j-fc, _|_ . . . 



where i > 1 and hj = (&2,j; • • • , ^nu)"^ is the vertical coefficient vector from 
(12.91) . For 1 < i < 111, comparing (12.101) with (12. 8p gives the recursive relation 
that computes the coefficients bj 

A^^\t,p2, . . .,Pm)bj = -«(}J„i],Ai+i(^'P2, • • . ,p„i,r„j+i, . . .) + lower terms, 
where 

'^{n' ,n],j [,(^n'+l,ji ■ ■ ■ i (^n,j ) ■ 

Then we substitute the formula for bj into the other Ui in (12.101) and get their 
updated Laurent r-series 



U^ = C^i{t)T-'^' + ■■■ + a^l_,{ty^-'-'^ 



.(2) 



,(2) 



+ ali^{t,p2, . . .,PmK' "' + ■■■ + all^_,{t,p2, . . .,Pm)r 



\2-l-ki 



'^ (^i,\2it, P2, ■ ■ ■ , Pni,1^ni+1, ■ ■ ■ ,fn2)'^ 

+ ■■■ + aS,_i(t, P2, • • • , Pm, rni+i, • • • , r^Jr^^"^"'"' 



+ ••• 

+ afi(^, P2, • • • , Pni, '^ni+i, . . . , r„Jr^^-^' + ■ 



(2.11) 



i > Hi, 



10 



where 

„(2) _ ^(1) „ / ^ 



(2) -J^)(f n n r \ 



^«J ""JJ 



l^V 2) • • • ) ' ni J 

+ lower terms, j > \i. 
Written in vector form, for j > Ai we have 

„(2) _.(!) _ nm(A(^))-^J^) p(l) _ ^"("i.»].Ai 

«(ni,n]j - «(ni,n],j ^ l^ ) «(l,ni],i' -« " ^(^ ) 



'•i=Pi 



,(2) 



For j = \i, I > 1, and n/_i < s < rii, the formula for a^ , . implies that 
the updated resonance vectors of resonance A^ form the matrix 

d{rni_i+i, ■ ■ ■ ,rni) 

a (1) a (1) 



^(^"i-i+l) • • • , ^nj C^(''n,_i + 1, • • • , rni) 



Therefore the updated resonance matrix 



is obtained as part of the row operation on R^^^ that uses the invertible 
matrix A^^' to eliminate B^^' 

I 0\ ^(1) fAW CW 



In particular, the invertibility of R^^^ implies the invertibility of R^'^\ 

The process continues and follows the same pattern. After introducing 
p„,_2+i, . . . , Pm^i, we get an updated Laurent r-series for the remaining Ui, 
i > ni_i. The coefficients of the series are functions of t,p2, . . . , Pm^i and 

11 



^n;_i+i) ■ ■ ■ T^n- Moreovei, we also have the resonance matrix R^^^^\ which is 
an (n — ri/_i) x (ri — ri/_i) invertible matrix. By rearranging the remaining 
Ui, we may assume that the first nii x mi matrix A^''~^^ in R^''~^^ is invert- 
ible. Then we introduce new variables by truncating the Laurent r-series of 
Uni_j^+i, . . . ,Uni at the resonance A^ and replacing the resonance parameters 
with new variables Pn^^-^+i, ■ ■ ■ , Pm- 



+ ■■• (2.12) 

The invertibility of A^''^ then implies that we can find the power r-series for 
the new variables 

pi = ri + I3i^i{t, P2, . . . , Pn,_i, ?^n,„i+i, ...,rn,)T^ 

+ /3i,A,+i-A,(t, P2, • • • , p„,_i, r„,_i+i, . . . , r„^+jT^'+'"^' + ■ ■ ■ 

+ ••• (2.13) 

+ /3i,A.-A,(t,P2,...,p„,_i,r„,_^+i,...,r„Jr^^"^' H , rii^i < i < rii, 

and the similar power r-series for the resonance parameters 

n = Pi + &i,i(t,p2,...,p„JrH 

+ &i,A,+i-A,(t, P2, • • • , p„,, r„,+i, . . . , r„Jr^2-Ai ^ 

+ ■■• (2.14) 

+ bi^Xs-\i{t, P2, ■ ■ ■ , Pni,rnt + 1, ■ ■ ■ ,rnjT^'~^' H , «;_! < i < Ui. 

Then we substitute (12.141) into the Laurent r-series of the remaining Ui, 
i > rii. We get the updated Laurent r-series, in which the resonance param- 
eters r„j_j+i, . . . ,rnj are replaced by the new variables Puf^^+i, ■ ■ ■ , Pnr The 
updated resonance matrix i?^') is obtained from the row operation on R^''~^^ 
that eliminates the terms below A'-'"^-'. 

We get the whole change of variable after exhausting all resonances. 



12 



2.3 Regularity 

Now we argue that the Laurent series in the balance (12. ip are converted to 
power series for the new variables r, p2, ■ ■ ■ , Pn- 

We have the power series (12.41) for r, with analytic functions of t, r2, . . . , r„ 
as coefficients. By taking the Taylor expansions in (f — fg) of the coefficients, 
we get the power series 

r = ai(to)(t - to) + ■ ■ ■ + ax,-i{to){t - to)^^"^ 

+ aAi(to,r-2,...,r„J(t-to)^' H 

+ ■■■ (2.15) 

+ ax^{to,r2, . . . ,rnj{t - to)^' H , 

with analytic functions of to, 1^2, ■ ■ ■ , ^n as coefficients. We note that ai(to) = 
i_ 

/3i(to) = cn,o(to) ^^ is nonzero. 

We have the power series (12.131) for pj. If n;_i < i < rii, then the co- 
efficients are analytic functions of t, p2, ■ ■ ■ , Pm^-^, fn,_i+ii • • • , ^n- We may 
substitute the power series (I2.15P for r and successively substitute the power 
series for pi corresponding to smaller resonances to the power series for pi 
corresponding to bigger resonances. Then we get updated power series for pi, 
with analytic functions of t, r2, . . . ,r„ as coefficients. Finally, we may take 
the Taylor expansions in (t — to) of these coefficients and get the power series 

Pi = Ti + ai,i(to,r2,...,r„J(t-to) H 

+ «i,A,+i-A, (to, r2, . . . , rn,^,){t - toy'+'-^' + ■■■ 

+ ■■■ (2.16) 

+ ai,\s-Xiito,r2, . . . ,rnJit - to)^"'^' H , ra,_i < i < ni, 

We conclude that the Laurent series for ui, . . . , m„ are converted to power 
series for r, p2, . . . , Pn- Moreover, the new variables satisfy the initial condi- 
tions 

T(to) = 0, T'(to) = ai,o(to)"^ ^ 0, p(to) = r,. 

Next we argue that the new system of differential equations for r, ^2, • • • , Pn 
is regular. The formulae (12. 3 p and (I2.12p for the change of variable imply 
that the new system is 

^' = fl'l(t,r,/72,---,Pn), p'i = gi{t,T,p2,...,pn), 

13 



with the right side 

9i = liit, r, P2, ■ ■ ■ , Pa) + 0i,l(t, P2, ■ ■ ■ , Pn)r~^ H h 0i,JV,(t, P2, • • • , Pn)T"^% 

(2.17) 
where 7^ and (^jj- are analytic in their variables. Since the Laurent series 
(12.11) is a formal solution of the original system, the power series (I2.15P and 
(I2.16P is a formal solution of the new system. Substituting the power series 
into the z-th equation, we find that the left side is a power series, while the 
right side is a Laurent series, with a^ '0i,Ar,(to, r2, . . . ,r„)(t — to)"^' as the 
lowest order term. Therefore we conclude that 

4'i,N,ito,r2,...,rn) = 0. 

Since this holds for all to and resonance parameters r2, . . . , r^i, and ^j^at. is 
analytic, we conclude that (pi^Ni is constantly zero. This completes the proof 
of the regularity of the new system of equations. 

2.4 Triangular Change of Variable 

Our change of variable is not quite triangular. To get a triangular change 
of variable, we keep the first variable r and modify the construction of new 
variables by introducing 

+ ■■■ (2.18) 

+ PiT^'~''\ ni-i<i<ni, 

instead of fl2:T2D . 

The variables pi are related to pi by 

~ _ (0 /+ N 

Pi — ^J,A,l'^'P2, • • • iPn,_i,Pn,_i + l, ■ ■ ■ i Prii ) 

+ cS,(^'^'P2,---,Pn,_i), rii_i <i<ni, 

where q { is analytic in its variables. The invertibility of A^''^ implies that 
the relation is invertible, and the inverse is also analytic. The variables pi 
satisfy the initial condition 

14 



Since A^''^ are invertible, the initial values can be any number. 

Finally, we remark that the proof for the regularity of the new system 
can be easily applied to general triangular change of variable 

ui = r~^ 

Ui = ai(t, T,p2,..., Pi-i) + bi{t, r, p2, • • • , Pi-i)pi, 1 <i <n. 

The key here is that, because fi are polynomial in Ui, the right side of the 
new system is of the form fl2.17p . Then the power series solutions for pi with 
all the possible numbers as the initial values imply that Qi have to be also 
analytic in r. 

The remark on the general triangular change of variable also shows that 
fi do not have to be polynomials. All we need is that the triangular change of 
variable converts the right side to be meromorphic in r. This often happens, 
for example, when /j are certain rational functions. 

3 Hamiltonian Structure 



Consider a Hamiltonian system 






, dH 


P' = 


dH 

dq' 


where 






q= (gi,...,g„). 


p = 


- (Pi,-- 



,Pn), 

and H = H{t,q,p) is analytic in t and polynomial in q and p. Consider a 
balance 

Qi = ai^oito)it - toy''' + ■■■ + aij{to, rg, . . . , r„J(t - to)^~'' + ■ ■ ■ , 
Pi = bi,oito)it - to)"'' + ■ ■ ■ + hj{to,r2, . . . ,r„J(t - to)^'~^' + ■ ■ • , 

where the coefficients for (t — toY'^' and (t — to)-'"'' depend only on the 
resonance parameters r2, . . . , r„; with resonance < j. 

The following is the exact condition for Theorem B. The definition is 
justified near the end of Section HI 

Definition. A balance of the Hamiltonian system is Hamiltonian principal, 
if the resonance vectors form a simplectic basis of R^", and there is d, such 
that ki + li = d — 1, and Xi + pi = d — 1 for the resonances A; and pi of 
simplectically conjugate resonance vectors. 
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The symplectic property depends on the order of vectors in the basis, 
and this order is not the same as the order of the column vectors for the 
resonance matrix. Consider the usual resonance matrix 



ix — [Ro Hi 



Rm 



where the columns of Ri have resonance A^. Since A; is increasing in /, the 
condition A^ + /i/ = d — 1 implies that a column vector in Ri is simplectically 
conjugate to a column vector in Rg-i- This is incompatible with our usual 
order in a simplectic basis vi, . . . , f2„, where Vi is simplectically conjugate to 
Vn+i- So we need to reverse the last n columns of R and get the matrix 



S = R 



In O 

O T„ 



T 



/O 




l\ 

1 



Vi ■■■ 0/ 
Then the definition requires that S" is a symplectic matrix 



iJ Jo — tVj fJ 



O In 
-In O 



Using the notation in the definition, we may denote the resonances by 
(s = 2s + l) 

-1 = Ao < Ai < ■ ■ ■ < As < yUs < ■ ■ ■ < ^1 < ^0 = c?, /"/ = Ms+l-l- 

Since the resonance vectors form a simplectic basis, the columns of Ri and 
R2S+1-1 are simplectically conjugate and therefore the two matrices have the 
same number of columns. In case Xg = /Xg, the corresponding resonance 
vectors actually form the matrix [Rg Rg+i), where the columns of Rg and 
Rg+i are simplectically conjugate and therefore the two matrices have the 
same number of columns. 

3.1 Simplectic Resonance Variable 

The proof of Theorem A may be adapted to prove Theorem B. We only need 
to be more careful in keeping track of the simplectic structure. 

Note that if the indicial normalization is assigned to gi, then the last new 
variable, which corresponds to the largest resonance A2S+1 = fJ-o = d, should 
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be introduced for pi. Therefore if we follow the construction in Section HI 
we should rearrange the variables in the order gi, . . . , g„, p^, ■ ■ ■ ,Pi. The 
resonance matrix corresponding to this order is obtained by reversing the 
second half of the rows of R 

D'= f^n 0\ (In 0\ (In O 

\0 Tn)^ \0 TnJ^'yO Tn 

The resonance matrix is updated by row operations. This is equivalent 
to an LU decomposition for R' , where L is lower triangular and U is upper 
triangular. In Section |2l the triangular shapes are actually blockwise. How- 
ever, we get LU decomposition only up to exchanging rows of the resonance 
matrix. For a Hamiltonian system, we need to exchange rows of R' so that 
the simplectic structure is still preserved. 

We divide the simplectic matrix into four n x n blocks 

The columns of j „ I span a Lagrangian of the simplectic space M?"'. It is 

well known in simplectic linear algebra that there is a subset / C {1, . . . , 2n}, 
such that 

1. The projection ttj: M^" — > M^ sends the columns of j „ J to a basis of 

w. 

2. For any I < i < n, one oi i,n + i is in J, and one is not. 

Now for any 1 < i < n, if n + z G /, then we exchange the z-th row and the 
{n + i)-th row of S, and ii n + i G /, then we do nothing. After this exchange, 
S is still a simplectic matrix, with A invertible. Correspondingly, this means 
that Qi and pi are exchanged. But the system of differential equations is still 
Hamiltonian under such exchange. 

So without loss of generality, we may assume that A is invertible. By 
further exchanging rows of A, we have A = LU. Exchanging rows of A 
means exchanging qi with Qj. This should be balanced by exchanging pi with 
Pj in order to keep the system Hamiltonian. Therefore if the i-th row and the 
j-th row of A are exchanged, then the i-th row and the j'-th row of S should 
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be exchanged, and the (n + i)-th row and the {n + j)-th row of S should also 
be exchanged. Such exchange keeps S to be simplectic. 

After all the "canonical exchanges", which keep 5* simplectic and the 
system Hamiltonian, we have a decomposition A = LU that fits the con- 
struction in Section [2l Then the fact that S" is a simplectic matrix implies 
the LU decomposition of the resonance matrix 

r., _( L O \(U L-^BTn 

\T^CA-^L t^{l^)-'tJ [o uu^y'T„ 

The construction of the change of variable in Section [2] is guided by this 
LU decomposition. The triangular version (12. 3p and (12.181) of the change of 
variable is 



qn = anfiqi^" + --- + an,x,-iqt ^ ''" + ^n^l ^ ''" 



Pn = Pnfiql'" + ■■■ + /3n,M.-l-;„?T ' + Pn^T 



(3.1) 



P2 = /32,ogr'^ + ■ ■ ■ + l32,,,-i-l,q^'~'-'' + p2q^'~'', 

Pi = /3i,ogT'^ + ■ ■ ■ + /3i,Mo-i-^i9T"'"'^ - K'Piqi'~'\ 



where aij is an analytic function of t,q2, . ■ ■ ,qi-i, and (3ij is an analytic 
function of t,q2, ■ ■ ■ , qn,Pn, ■ ■ ■ ,Pj+i- Note that a numerical factor — A;f ^ is 
added in the last line. 

3.2 Canonical Change of Variable 

In this part, we argue that the change of variable is canonical 

dqi A dpi + ■ ■ ■ + dqn A dpn = dqi A dpi + ■ ■ ■ + dijn A dpn- 
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From the shape of the triangular change of variable fl3.ip . we have 

dqi = -kiq^^^''^dqi = -kiql^'^^dqi, 
dpi = — A;f ^gf"^ ^dqi + linear combination of qf^°~ ^ , 
dqi = qi~ 'dqi + linear combination of gf '~ % 
dpi = g^'~ "dpi + linear combination of gf^'~ \ 

By {Xi - ki) + {fii - k) = {Xi + fii)-{ki + k) = {d-l)-{d-l) = 0, we have 

dqi A dpi + ■ ■ ■ + dqn A dpn = dqi A dpi + ■ ■ ■ + dqn A dpn 

+ q^^ioi + ■■■ + qi^UJN, (3.2) 

where Ui is a 2-form in gi,...,g„, pi,...,pn, with analytic functions of 
t,q2, . ■ ■ , (jn, Pi, ■ ■ ■ ,Pn as coefficients. 

Think of the Laurent series in the Hamiltonian principal balance as a 
transform between gi, . . . , g„, pi, . . . ,pn and to,r2, . . . , r2n- We substitute 
the transform into the left side of (13.21) and get a 2-form in to, ^2, . . . , r2„, 
with functions of t, to, ''^2, • • • , ?'2n as coefficients. Since the symplectic form is 
invariant under the Hamiltonian flow, the coefficients are actually indepen- 
dent of t. 

On the other hand, the Hamiltonian principal balance is converted to 
power series for the new variables qi, ... ,qn, pi, ■■■ ,pn- Think of these power 
series as a transform between gi, . . . , g„, pi, . . . ,p„ and to,r2, . . . , r2„. We 
substitute the transform into the right side of (13. 2p and get a 2-form in 
to, r2, . . . , r2n, with functions of t, tg, r2, . . . , r2„ as coefficients. Since this 2- 
form is the same as the substitution of the Laurent series into the left side, 
the coefficients of this 2-form are independent of t. 

Let 

UN = dqi A p + (f), 

where p and (j) are 1-form and 2-form in g2, . . . , g„, Pi, . . . ,pn, with analytic 
functions of t, g2, . . . , g„, pi, . . . ,pn as coefficients. After substituting the 
power series into un, we get 

ujn = dto A p + (p + {t - to)?/), 

where 

1. p and i/j are 1-form and 2-form in r2, . . . ,r2n, with analytic functions 
of t, to, r2, . . . , r2n as coefficients. 
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2. (j) is obtained as follows: Think of the initial data as a transform 
(r2, . . . , r2n) ^-> (92(^0), • • • , qn{to),Pi{to), ■ ■ ■ ,^n(^o)) from the resonance 
parameters to new variables, with to as a constant. Substituting the 
transform and t = to into (j) gives (j). 

The shape oiuj^f implies that, after substituting the power series for gi, . . . , g„, 
pi, . . . ,pn into the right of fl3.2p . we get a 2-form 

Here ai is the leading coefficient in f l2.15p for t = qi, and depends only on to 
and resonance parameters rj of resonance 0. Moreover, ■ ■ ■ consists of terms 
that are either of order > — A^ in (t — to) or have dto as a factor. Since we 
have already argued that the coefficients of the 2-form should not depend 
on t, and a^^ is independent of t, and (j) is independent of t and has no dto 
factor, we conclude that = 0. Since is obtained from (p by an invertible 
transform, this implies = 0. 
So we have 

U]\f = dqi A p = dtf) A p + {t — to)'0. 

Since there is no more 0, we now know that p is obtained from p in the 
similar way as being obtained from above, and then multiplied by —ai. 
Moreover, -0 has the same description as above. Therefore the right side of 
(13.21) now becomes 

a^^{t-to)-^dtoAp + --- , 

where ■ ■ ■ consists of terms that are of order > — A^ in (t — to). Since the 
coefficient of the 2-form should not depend on t, and ai and p are independent 
of t, we conclude that p = 0. Since p is obtained from p by an invertible 
transform, this implies p = 0. 

We proved that Ui\f = 0. Therefore there is no negative powers on the 
right of fl3.2j) . so that the change of variable f l3.ip is canonical. 

The argument was carried out for the triangular change of variable. It 
also works for the blockwise triangular change of variable, making use of the 
simplectic property of the resonance matrix. Alternatively, we may verify 
(again by the simplectic property of the resonance matrix) that the transform 
between the new variables in the blockwise triangular change of variable and 
the new variables in the blockwise triangular change of variable is canonical. 
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3.3 Change of the Hamiltonian Function 

In this part, we study the Hamiltonian function for the new Hamiltonian 
system. 

Suppose the original Hamiltonian system is autonomous. Since our change 
of variable is canonical, the new system for the new variables is also an au- 
tonomous Hamiltonian system, and the new Hamiltonian function is simply 
obtained by applying the change of variable to the original Hamiltonian func- 
tion. 

For an autonomous (not necessarily Hamiltonian) system, the coefficients 
in a balance do not depend on t. Then in the construction in Section |21 all the 
coefficients in the Laurent and power series are independent of t. Therefore 
the regularizing change of variable we get at the end is independent of t. 

For the autonomous Hamiltonian system, applying the triangular change 
of variable (13. ip to the original Hamiltonian function, we find the new Hamil- 
tonian function to be of the form 

H{q,p) = ho + hiq^^ H h HnQi^ , 

where ho is a polynomial of gi, . . . , g„, pi, . . . ,pn, and hi, ... .h^ are poly- 
nomials of q2-, ■ ■ ■ iQni Pi, ■ ■ ■ iPn- Since the Hamiltonian function is a first 
integral of the Hamiltonian system, substituting the power series solution for 
qi, . . . ,qn, pi, ... ,pn into H should give us an expression independent of t. 
The most singular term after the substitution is 

hN{q2{h)-: • • • , qn{h),Pi{h)-, ■ ■ ■ ,pn{to))ai^{t - to)"^- 

Since hj^ and ai are independent of t, we conclude that 

hN{q2{to), ■ ■ .,qnito),pi{to), . . . ,Pn(^o)) = 0. 

Since the initial value map (r2, . . . , r2„) ^-> (52(^0), • • • , qn{to) , pi{to) , . . . ,pn(io)) 
is invertible for triangular change of variable, we conclude that h^ is a con- 
stantly zero function. This proves that H{q,p) = ho is a. polynomial. 

Now we turn to non-autonomous Hamiltonian system. In this case, the 
function H{t,q,p) obtained by applying the time dependent change of vari- 
able to the original Hamiltonian function H{t,q,p) is no longer gives the 
Hamiltonian function of the new system. In fact, H may be singular in qi. 
In what follows, we argue that it is the "regular part" hQ{t, q,p) that becomes 
the Hamiltonian function of the new system. 
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We leave Hamiltonian systems for a moment and consider a general sys- 
tem of differential equations 

U' = f{t,u), U = {Ui,...,Un). 

For a change of variable u = ip{t, p), we have 



u 



dt ^ dp^' 



So the new system is 



p' = J-'f{t,^{t,p))-J-^^, J, = g. 



For the triangular change of variable (12.31) and fl2.18p {p = (r, p2, • • • , Pn) 
takes the place of p) , we have 



/-kiT' 



i-fci 



J. 



-<Ai-fc2l 



T 







\ 




\{- 



<, As Kn 



I r^<As-fc„| 



T 



^s <^n 



where {t^-'} means a (finite) linear combination of t^ with f < j, with 
analytic functions oi t, p2, ■ ■ ■ , Pn as coefficients. Moreover, /(t, ip(t, p)) and 

-— are also (finite) linear combinations of r-', with analytic functions of 

t, p2, . . . , Pn as coefficients. Then both J~^f{t, ip{t, p)) and J~^—- are linear 
combinations of t^ of the same type. 

Let J^^-^ = {9i,---,gnV- Then — = J^{gi, ■ ■ ■ ,gnV, and by the 
formulae (12. 3p and (I2.18p . we get 



i-fei 



= -k,T-'-'''g,, 



|^<Ai-fc2}_|^<Ai-A,-2|^^^^ 



<Ai-fc2^ 



rAi-A;2 



92, 



-\As /i^TT, 



|^<A,-fe„| ^ |^<A,-fe„|^^ ^ {r<^='^"}^2 + ■ ■ ■ + r 



<As-/i;nl 



^As fejT, 



22 



By induction, it is easy to see that gi = 0, g2 = {t^°}, • • • , gn = {'^^^}- 

_, dip 
Therefore we conclude that J,„ -:r- is a hnear combinations oi t\ j < 0, 

^ dt ^ 
with analytic functions oi t, p2, ■ ■ ■ , Pn as coefficients. Since the differential 

equation 

obtained after the triangular change of variable ( 12. 3p and (12.18^ is regular, we 
conclude that J^^-— is the sum of terms in J~^f{t,(p{t,p)) with r-', j < 0. 
In other words, the new equation is 



p'=[J-'f{t,^{t,p))] 



regular ) 



with the right side obtained by dropping terms with negative power of r. 

Back to the non-autonomous Hamiltonian system, since the triangular 
change of variable (13.11) is canonical, we have 



j: 



dH \ 
H 



% 



dq I 



( d_H_\ 

dp_ 
dH 

\~~d^J 



Here //(t, q,p) is obtained by applying the triangular change of variable (13. ip 
to H{t,q,p) and is therefore a linear combination of ql, with analytic func- 
tions oi t,q2, ■ ■ ■ ,qn, Pi, ■ ■ ■ ,Pn as coefficients. We separate the parts of H 
with (fi, j >0 and the rest part with ql, j < 0, 



H = [-f^] regular + [H 



singular • 



Then it is easy to see that 

dH _ 
dp 



'dH 

dp 


+ 

regular 


'dH 

dp 



singular 



with 



dH 

dp 



regular 



(-/[il J regular 

dp 



dH 

dp 



d[H] 



singular 



singular 



dp 
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pjTT 

and the same happens to — -r- Thus we conclude that the new system of 



equations is 



q = 77^ — , p — 



dp dq 

This shows that [i/] regular is the Hamiltonian function of the new system. 

4 A Mathematical Foundation for the Painleve 
Test 

Let fi{t,Ui, . . . , Un) be analytic in t and polynomial in mi, . . . , m„. We consider 
a system of ordinary differential equations 

u[ = fi{t,Ui,...,Un), l<i<n. (4.1) 

The Painleve test attempts to find all Laurent series solutions (or balances) 
with movable singularity to 

Ui = dit - toy''^ + ai,i(t - to)^"''" + ■ • ■ + aij{t - toY'^'^ + ■■■ . (4.2) 

Here we use q instead of aj^o to highlight the different role from the subse- 
quent coefficients. 

4.1 Dominant Balance 

The first step in the Painleve test is to determine the leading exponents 
ki and the leading coefficients q of potential balances. There are usually 
several possible combinations of leading exponents and leading coefficients, 
which give several possible balances. 

Definition. The leading exponents /ci, . . . , A;„ of a balance is Fuchsian, if the 
A;*-weighted degree of /« is < /cj + 1. 

The k^-weighted degree of a polynomial in mi, . . . , w„ is obtained by taking 
the degree of Ui to be k^. 

If all the leading coefficients q 7^ 0, then the choice of leading exponents 
is natural. For the natural leading exponents A;*, denote the dominant part 
of/. 

Dj(t, wi, . . . , Un) = y ^ terms in fi with highest A;*-weighted degree. 
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The following is a simple criterion for the natural leading exponents to be 
Fuchsian. 

Proposition 1. // Dj(to,Ci, . . . , c„) 7^ and kiCi ^ for each i, then the 
natural leading exponents is Fuchsian. 

Let di be the /c*-weighted degree of /j. Then after substituting the balance 
f l4.2p . the lowest order term in u'^ is —kiCi{t — tg)^^'^^, and the lowest order 
term in /j(t, Ui, . . . , Un) is Di{to, Ci, . . . , c„)(t — to) *• Under the assumption 
of the proposition, both have nonzero coefficients. Therefore they are the 
"real" lowest order terms, and we conclude that 

di ^ ki -\- i, Ui(tQ, Ci , . . . , Cn) = fcjCj- 

The first equality shows that the Fuchsian condition is satisfied. 

If the natural exponents is not Fuchsian, then we need to choose some- 
what unnatural leading exponents in order to satisfy the Fuchsian condition. 
For example, the Gelfand-Dikii hierarchy with 2 degrees of freedom is a 
Hamiltonian system with 

H = -qivl - 2pip2 + 3qfq2 - qt - ql- 

One of the principal balances of the system is (we abbreviate (t — to) as t 
because the system is autonomous) 

^1 = t'^ + 3^2 - -rlt" - -rst^ - —rlt" - -r^r-^t'' - -r4^ + . . . , 

92 = r^t-^ - \rl - r;t - K^t^ + (-^4 + ^rA t" + ■ ■ ■ , 

_, 2 _. 1 4.5 2 /'22 4 11 \ . 

Pi = -t + -r2t -^ + -ra rU r2r3t^ + — r^ r^ f H , 

Fi 3 6 27 ^ 6 V81 3 7 

P2 = t"^ + -rlt + rst'^ + —rlt^ + -r2rzt'^ + Tit'' + .... 
6 lib 

The natural leading exponents 2, 2, 5, 3 is not Fuchsian. By taking the leading 
exponent of 

52 = Ot-^ + 0t~3 + r2t"2 - %-l - rst - ^rlt^ + (-^4 + ^^4] t^ + • • • , 
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to be 4, we get leading exponents 2, 4, 5, 3 satisfying the Fuchsian condition. 
The corresponding leading coefficients are 1, 0, —1, 1. The Henon-Heiles sys- 
tem in [12] is another example. 

When the Fuchsian condition is satisfied, we denote the dominant part of 

fi 

ff = N^ terms in fi with ^^.-weighted degree ki + 1. 

Under the condition of Proposition [T], we have ff = Di. By substituting 
the balance (14. 2 p into the system (14.11) and comparing the coefficients of 

(t — to)"'''~\ we get 

f[^{to,Ci,...,Cn) = -kid. (4.3) 

This is the (formal) dominant balance for the system fl4.ll) at the (movable) 
singularity to- 



4.2 Kowalevskian Matrix 

Given the leading exponents and the leading coefficients, the next step of 
the Painleve test is to find the subsequent coefficients in the balance. For 
j > 0, by substituting the balance (14. 2 p into the differential equations and 
comparing the coefficients of (t — to)-'"'''"^, we get 



(j -ki)a. 



dfF 



*j 



(^0) Ci, 



df, 



D 



' ''"^''''^' ^'"^ du. 



dui 
+ terms involving c*, a^K^i 



\''0) '-'1) ■ ■ ■ ) (^n)0'n,j 



''^^j — 1 • 



Define the Kowalevskian matrix 



dui 



K 



ki 



dJl 

dui 



dfl 
du2 

ir^ + ^2 

OUo 



D 



\ dui 



duo 



D 



dUF. 



d{ui,. 



-(to,ci, 



dUr, 



,Ur 



dUr. 



D 



dfl 

dUn 

(ki 



+ fen 



■■,Cn) + 



\ 
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Then we have a recursive relation 

^"^ — terms involving c*, a* i, 



{K-jI){a,„ 



, a 



n,]) 



I ^*,7 — 1 • 



(4.4) 



A formal balance is the Laurent series solution with leading coefficients 
obtained by solving ( I4.3p and with subsequent coefficients obtained by solving 



We emphasize that the Fuchsian condition is needed for the definition of 
the Kowalevskian matrix. Ercolani and Siggia |5] made a similar observa- 
tion regarding the balance for the Gelfand-Dikii hierarchy. They suggested 
that one should not blindly proceed to use the linearized equation as the 
constraint on the resonance vectors and specifically pointed out that the 
leading exponent for q2 should be changed from 2 to 4 in order to define the 
Kowalevskian matrix, which is 



K 



/ 2 





-2^ 


-2 4 


-2 


-2 


12 -6 


5 


2 


\-Q 2 





3/ 



The solutions of fl4.3p and fl4.4p may not be unique. This leads to free 
parameters in the balance. Since (14. 4p is a linear equation, we can choose 
the free parameters such that if a parameter does not appear in the leading 
coefficients, then the parameter first appears linearly. Although this property 
is not needed for our main theorems, it is a consequence of our set up. 

There are three kinds of free parameters. The first kind is the free location 
to of the singularity. The derivative of the balance (14. 2 p in to is 



dtf) 



dc' 

[1 ~ ki]aii + -—^ 
' ' dto 



{t-to)-'^ + 



Therefore the variation of the solution due to the parameter to is character- 
ized by the basic resonance vector {—kiCi, . . . , —knCn)- 
Take the derivative of the system (14. ip in t, we get 

,/ 9fi .dfi, dfi , 

Substituting in the balance (14. 2 p and comparing the coefficients of (t — 

to)'^''"^, we get 



{-ki){-ki-l)ci 



-(to, Ci, . . . , Cn)\~ki)Ci + - ■ ■ + ^ (to, Ci, 



dui 



dur 



I ^n ) \ '^n ) Cn • 
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The equalities can be rewritten as {K + J)(— fciCi, . . . , —knCn) = 0. 

Proposition 2. // the leading coefficients of a balance ( 14. 2 p is Fuchsian, 
then the basic resonance vector (— fciCi, . . . , —knCn) is an eigenvector of the 
Kowalevskian matrix with eigenvalue —1. 

The second kind of free parameters appear in the leading coefficients, and 
parameterize the subvariety given by f l4.3p . The derivative of the equation 
( 14. 3 p with respect to one such free parameter r is 

^^^dc, _df['dc, ^ ^ dfPdc^ 

Or dui dr dun dr 

The equalities can be rewritten as K ' ' " — = 0. 

or 

Proposition 3. // the leading coefficients of a balance (14. 2 p is Fuchsian, 
then the tangent vectors to the subvariety of allowable leading coefficients are 
the eigenvectors of the Kowalevskian matrix with eigenvalue 0. 

The third kind of free parameters appear in the subsequent coefficients. 
These are obtained from solving the recursive relation (14. 4p . One has to 
worry about the compatibility between K — jl and the right side of the 
recursive relation, which affects the existence of solutions. Leaving aside 
the existence issue for a moment, we may conclude the following from the 
recursive relation. 

Proposition 4. // the leading coefficients of a balance (14. 2 p is Fuchsian, 
then for j > 0, the j-th coefficient vectors form an affine space parallel to the 
eigenspace of K with eigenvalue j. 

If j is not an eigenvalue of K, then the j-th coefficients are uniquely 
determined by (14. 4p . If j is indeed an eigenvalue, then the j-th coefficients 
form a vector 

(aij, . . . , a„j)'^ = ttj + riRi -\ h TmRm, 

where 

1. cij depends only on to, c*, a,,^i, . . . , a*j_i. 

2. Ri, . . . , Rm form a basis of the eigenspace of K with eigenvalue j and 
therefore depend only on to and c*. 

28 



3. ri, . . . , r^ can be any numbers. 

Definition. A Fuchsian leading exponents is principal, if for all solutions of 
f l4.3p . the following are satisfied. 

1. The Kowalevskian matrix K is diagonalizable, with the eigenspace of 
— 1 being of dimension 1, and all other eigenvalues being non-negative 
integers. 

2. The recursive relation (14. 4p is always consistent. 

Given principal Fuchsian leading exponents, the corresponding principal 
balance (14. 2 p is obtained by solving equations (14. 3 p and (14. 4p . And the 
Painleve test is passed. 

The first condition implies that the eigenvalues of K and their multiplic- 
ities are constants. By applying the implicit function theorem to (14. 3p . we 
also know that the collection of leading coefficients form a submanifold with 
the eigenspace ker K of eigenvalue as the tangent space. Therefore the 
definition means exactly that the total number of free parameters, including 
^0, is n. Moreover, the eigenvectors of K form the resonance matrix. 

For the balance of the Gelfand-Dikii hierarchy, the Kowalevskian matrix 
has eigenvalues —1,2,5,8, with the corresponding resonance matrix 

/ 2 1-4 -2 \ 

3-69 

-5 2 1 -22 

\ 3 6 6 / 



R 



4.3 Painleve Test for Hamiltonian System 

Consider a Hamiltonian system given by a Hamiltonian function H{t,q,p) 
that is analytic in t and polynomial in q and p. For a balance, we denote the 
leading exponents of Qi and Pi by ki and k. 

Definition 5. A Hamiltonian system is almost weighted homogeneous rela- 
tive to the leading exponents ki, . . . ,kn, h, ■ ■ ■ Jn ii ki + k = d — 1, where d 
is the leading exponents weighted degree of H. 

The almost weighted homogeneous condition implies that the weighted 

degree of — — is < d — ki = li + 1 and the similar inequality for — — . There- 
opi dqi 

fore the Fuchsian condition is satisfied. Let H^ be the sum of those terms 
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in H with highest weighted degree and apply the same notation to other 



polynomial functions. Then 
is 



K 



dqidpi 



dpi 



OH 

dpi 

dpi 

dpidqij y 



■, and the Kowalevskian matrix 
\ 



h 



•••/ 



O I 
-I o 



JU + T. 



dqf 



dpidQi 



(h 



+ 



\dqidpi dpi J 



h 



■J 



where "H is the Hessian of H^ . The special shape of the Kowalevskian matrix 
implies certain symplectic structures among the resonance vectors. 

Proposition 6. Suppose the Hamiltonian function is almost weighted ho- 
mogeneous relative to the leading exponents of a balance. Suppose v and w 
are eigenvectors of the Kowalevskian matrix with eigenvalues A and jj,. If 
A + /i 7^ d — 1, then {v, Jw) = 0. 

The condition ki + k = d — 1 means exactly TJ + JT = {d — 1) J. Then 
the lemma follows from 

(A + /i)(v, Jw) = {Kv, Jw) + {v, JKw) 

= -{j{jn + T)v, w) + {v, j{jn + T)w) 

= {{H - JT)v, w) + {v, {-n + JT)w) 
= {v, TJw) + {v, JTw) 
= {d — l){v, Jw). 

The proposition implies that, by rescaling the vectors if necessary, we can 
find a simplectic basis of resonance vectors. For example, by rescaling and 
exchanging the last two columns of the resonance matrix for the balance of 
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the Gelfand-Dikii hierarchy, we get a symplectic matrix 



R 



The last two columns are exchanged because the proposition tells us that in 
the original resonance matrix, up to rescaling, the first and fourth columns 
are simplectically dual to each other, and the second and third columns are 
simplectically dual to each other. 



A 


0\ 




(^ i 


2 

81 


-f\ 


1 







1 


1 
9 


^3 





\ ^ 




-5 1 


22 
81o 


1 


\Q 


-k V 




V 3 


2 

~27 


3 / 
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